Two Graviton Production at e + e and Hadron Hadron Colliders 

in the Randall- Sundrum Model 



Pankaj Jain and Sukanta Panda 

Physics Department 
IIT, Kanpur, India 208016 

(N : 
c : 

| c3 . Abstract We compute the pair production cross section of two Kaluza Klein modes in the 

Randall- Sundrum model at e + e~ and hadron hadron colliders. These processes are interesting 

^ | because they get dominant contribution from the graviton interaction at next to leading order. 

' Hence they provide a nontrivial test of the low scale gravity models. All the Feynman rules at 

^ ■ next to leading order are also presented. These rules may be useful for many phenomenological 

CN . applications including the computation of higher order loop corrections. 
(N ■ 
(N ■ 

O '. 1 Introduction 

^ The proposed existence of large extra dimensions |2| might provide a solution to the hierarchy 
*p ■ problem since it can considerably lower the scale of quantum gravity. In this model the standard 
± , model fields are assumed to be confined on a 3 + 1 dimensional hypersurface or a brane in a 
D \ higher dimensional space-time [B] . In this low scale gravity model, which is generally refered to 
^ \ as the ADD model, gravity may become strong even at TeV scale which is far below the Planck 
# ^ ■ scale 10 19 GeV. The two scales are related via gauss law given by = M^ n V n , where V n 
^ ■ is the volume of n extra dimensional space. If we take Mp ~ TeV then the size of the extra 
^ , dimension lies in the range 1 mm - 10 fm for n = 2 - 6. However by doing so we introduce a 
new hierarchy between the compactification scale and the fundamental scale in the theory. An 
alternative low scale gravity model (the RS model) [3] involves the existence of a single extra 
dimension in a warped geometry. In this case the scale of quantum gravity can be as small as 
1 TeV without requiring the existence of very large extra dimensions. These low scale gravity 
models are reviewed in Ref . El • The leading order Feynman rules for the ADD model are 
given in Lykken et al (7j and Giudice et al [H]- 

Phenomenologically the RS model jU] is very different from ADD model because of the 
presence of very massive KK modes which can produce observable signatures in future collider 
experiments. In the case of ADD model we instead have large number of closely spaced KK 
modes with the mass the lightest mode governed by the size of the large extra dimension. 
These modes cannot be observed directly due to their small coupling to matter and hence the 
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production of these modes gives rise to missing energy signatures. However in the RS model 
the first few modes produced can decay into observable particles. The center of mass energy 
required to produce these resonances are well within the energies available in the next generation 
of colliders such as the LHC Non-observation of these processes in experiments constrain the 
parameters of this theory. Several phenomenological studies have been carried out to study 
these resonances in e + e~ colliders (TU1 HH H21 Ed] and in hadron colliders |T%1 Hoi ITo] IT71 ITS]. 

In the present paper we study the process of two graviton production in the RS model in 
e + e~ and hadron hadron colliders. This process gets contribution at next to leading order in 
the gravitational coupling. Hence we derive the required Feynman rules at this order. This 
process may not be very important in case of ADD model because these modes just escape 
through the detector undetected. However our calculation can be easily generalized for this 
case also. 

The plan of the rest of the paper is as follows. In the second section we briefly review 
the derivation of lowest order interaction Lagrangian. We then extend the calculation to next 
order in coupling and find the corresponding interaction Lagrangian. Third section describes 
all processes related to emission of two lowest massive graviton modes. Here we also present 
our results for the cross section of two graviton emission at e + e~ and hadron hadron colliders. 
Section four gives our conclusions. In Appendix A we present all the Feynman rules involving 
vertex of matter fields with two gravitons at next to leading order in coupling. In Appendix 
B, we derive the Feynman rule for three graviton vertex which also contributes to the process 
under consideration. 



2 Interaction Lagrangian 

The RS model assumes a five dimensional non-factorizable geometry where two 3 — branes, 
with opposite tensions, reside at S 1 /Z2 orbifold fixed points. In this framework solution to 
5-dimensional Einstein equations is given by 

ds 2 = e-^rj^dx^dx" + r 2 c d<f) 2 (1) 

where is the angular coordinate parameterizing the extra dimension, cr(0) = kr c \(p\ with r c 
being the compactification radius of the extra dimension, and < \<p\ < n. Starting from the 
5-D action, after integrating over extra coordinate 0, we obtain the reduced effective 4-D planck 
scale given by the relation 

M 3 

M 2 = — (l- e - 2kr ^) . (2) 
k 

Setting visible brane at <p = it, where our standard model fields live, it is found that any 
mass parameter m on the visible 3-brane in the fundamental higher- dimensional theory will 
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correspond to a physical mass 



m 



e krc7T m . 



7r due to exponential factor 



Hence TeV mass scales can be generated on the 3-brane at 
present in the metric if we assume kr c « 12. 

A small perturbation of the metric, such that e~ 2a ^r]^ L , — > e~ 2a ^ (r/^ + h^), can be 
expanded as, 



n=0 



(3) 



After solving the linearised equation of motion of graviton field and working in the gauge 

d»h^) = h^ n) = 0, 

one obtains the solution for (</>)> 



X (n) (0) 



+ a n Y 2 



(Mn 

V k 



(4) 



The exact expression for a n is given in Ref. jSj. The masses of the KK modes, M n , are given 
by EH] 

M n = x n fce- fcrcW (5) 

where X^S £1X6 the solution of the equation J\(x n ) = 0. We also define the mass parameter 
itlq = M n /x n . As we see from Eq. masses of the graviton KK excitations are dependent on 
the roots of J\ and are not equally spaced. In the following, we formulate the interaction of 
physical KK modes to the Standard Model fields which reside on the brane at = ir. 

To start with, let us first consider the minimal gravitational coupling of the general scalar 
S, vector V and fermion F, 

J d*x^£(g,S,V,F) (6) 
where the induced metric g on the brane can be decomposed as, 

By expanding the interaction Lagrangian upto 0(k 2 ), we get 



d 4 x 



gC(g) 



d 4 x 



7T 



AC(g)\ 



Wh^(x,(j)) 



£(g)\v=v-^T, u 

5C 

Sg^ 

S 2 C 
bg^ v bgP u 



9=V 



\9=V 



(7) 
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where we have used 



1 2 ' 1 '- 2 \K„hr 



g = l + -h + K 2 (-h 



Several coefficients that appears in Eq. (JJJ) are given by, 

A = l h2 ~ \W , 

= -hh» v - h^hK , 
2 A 

h = h* = h^ Vfll/ . 
Interaction Lagrangian upto 0(k) is given by 



(9) 



(10) 

(11) 
(12) 



h tiV (x,(f) = tt)T. 



1 



fiv 



h^(x,4> = ir)T. 



/IV 1 



where k 



M372 

terms of KK modes as, 



. With the help of Eq. |2] we can express the above interaction Lagrangian in 



I 1 oo 

M Pl ^ A. 2 - ^ 



(13) 



l i" n=l 



where TL,, the energy- momentum tensor of the matter field confined on the brane, is given by, 



5C 



1 ixv 



\9=V 



(14) 



and 



A w = M Pl e- kr ^ . 

Now we are interested in finding out the 0(k 2 ) Feynman rules for the gravitational interaction 
with matter confined on the D 3 brane embedded in the RS bulk. 

The k 2 part of the Lagrangian that resembles the gravitational interactions with matter 
involving two gravitons, is given by 
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E E h^ n X n (<P = *)hTx m ^ = tt) 

n=0 m=0 

oo oo 

(1/2) E E ^X n (0 = n)h^ m X m ^ = 7T) 



5C 



5gt MU 



n=0 m=0 



9=V 



9=V 



oo oo 



+ (1/8) E E h n X n (<f> = *)h m X m (<f> = tt) C(g)\ S =r, 

n=0 m=0 

oo oo 

(1/4) E E K^ n (<P = vr)^ m X m (0 = tt) £^1^ 



n=0 m=0 

oo oo 



+ (V2) E E ^' n X n (0 = tt) „ t 
„=o m =o d^d^ 



9=V 



(15) 



After substituting the expression of x n i n the above equation and using the relation given in 
Eq. 121 we finally get 

£ 2 = £°° + 2£ 2 0n + C n 2 m (16) 
where the £™ part of the lagrangian is given by 



prim 
t-2 



^ oo oo 

T2 E E 

tt n=l m=l 



5C 



5C 



tig 



g=v 



g=v 



-(l/2)/W' ? , 
v ' ' 8g^ 

+(l/8)/W» C{g)\ s=r) 
-(mK„h»°> m C(g)\ s = v 



Sg^Sg 



pa 



9=V 



(17) 



The £9,° P ar t °f the lagrangian is obtained by setting n = m = and by replacing the coefficient 
A 2 n by M 2 Pl in £™ m . Due to the large suppression factor 1/M|, ; the contribution of this lagrangian 
as well as that of £2°™ is negligible. By using the interaction Lagrangian, we can derive the 
feynman rules for two gravitons interacting with matter fields which are given in Appendix A 
of this paper. Also the feynman rules for three graviton vertex are presented in Appendix B. 
The triple graviton vertex has also been studied in Ref. 



3 Two Graviton Production 

We next consider the production of two gravitons at e + e~, pp and pp colliders. We first consider 
the process fermion anti-fermion annihilation into two gravitons. There are four feynman 
diagrams contributing to this which are shown in Fig. ([T|. The differential cross section for 



5 



(a) (b) (c) (d) 



Figure 1: All possible tree level Feynman diagrams which contribute to two graviton emission 
due to fermion anti-fermion annihilation. The solid and dotted lines represent the fermions and 
the gravitons respectively. 

e + e" — > GG is given in Eq. EH] in Appendix C. The cross section also gets contribution from 
the s-channel diagram which involves a graviton propagator. This gets contribution from all 
the Kaluza-Klein modes. By summing over all the contributions the propagator takes the form 

32ttc| 1 32^ 

{ ) ml ^s-Ml+iY n M n mi MXs) ' ^ 

where x s = ^ and Cq = k/M p \. By taking into account all decay channels, T n can be written 

as 

T n = c 2 m x 3 n A(M n ) (19) 

Here A is a complicated function of M n . We calculate the A s (x s ) numerically and it's behaviour 
is shown in Fig El 

In Fig. El we show the s dependence of the production cross section of the two lightest KK 
modes of mass M\ = 3.83m in e + e~ collisions. The resonance peaks which arise due to the 
s-channel graviton exchange (diagram |TJi) are clearly visible in a certain range of parameter 
space. The peaks tend to disappear for large values of the parameter Co = k/M v \. We expect 
that the first few peaks will be identifiable if the parameter cq is sufficiently small. For values 
of Co larger than approximately 0.6 even the first resonance may not be clearly identifiable |T3*] . 
The maximum value of the parameter mo that can be explored with this process is given by 
y/s < 2M\ = 7.66m . In Fig. HI we show the angular dependence of differential cross section for 
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Figure 2: The graviton propagator factor A s as a function of x s = yfs/rriQ. Here we have fixed 
the parameters m = 100 GeV and c = .01. 

pair production of the first graviton mode with mass Mi = 3.83mo where mo has been chosen 
to be 150 GeV. 

We see from fig. El that, as expected, the cross section rises very rapidly with s. The 
perturbation theory is applicable only if s << A^. In our calculations we have chosen the 
parameter space such that this condition is respected and perturbative predictions can be 
trusted. Perturbation theory breaks down for values of Co much larger than those chosen in fig. 
E] for the corresponding values of m . In this region it is not possible to compute this process 
reliably with our current understanding of quantum gravity. 

The hadron-hadron crosss section gets contributions both from quarks and gluons. In the 
case of quark anti-quark annihilation the contributing diagrams are same as listed in Fig. d 
The diagrams contributing to gluon gluon fusion are also the same as those given in Fig. ^ 
with the quark lines replaced by gluons. The results for the pair production of the lightest KK 
mode with mass Mi = 3.83mo is given in figures Hand El The center of mass energy yfs of 
the hadron hadron system is taken to be the LHC energy, i.e. 14 TeV. In fig. E]we show the W 2 
dependence of the cross section, where W 2 = ts = X1X2S is the invariant mass of the final state 
and x\ and xi are the longitude momentum fractions of the two initial state partons. We see 
from this figure that, as expected, LHC can probe a considerable range of parameters of the RS 
model through this process. In figs. E]and[7|we show the cos# and the rapidity y = \ ^og{x\/x2) 
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Figure 3: The two graviton production cross section in e + e~ collisions as a function of the 
center of mass energy squared s (TeV 2 ) for several different values of the parameters (co, Too), 
where mo is given in GeV. 

dependence of the cross section respectively. Here 9 is the center of mass scattering angle in 
the parton subsystem. In Figs. |H1 El and EH we show the corresponding results for the proton 
anti-proton system at center of mass energy y/s = 2 TeV, relevant for the Tevatron. 

Each of the two gravitons produced in the final state eventually decay into either two jets, 
a lepton pair or a vector boson pair. Therefore one possible experimental signature of this 
process involves a four jet event such that the invariant mass of the two pairs of jets is equal 
to the mass of the graviton. There will be considerable QCD background which may limit the 
parameter space that can be explored with this process. We postpone the detailed study of 
experimental signatures to future research. 



4 Conclusions 

We have computed the pair production cross section of Kaluza Klein modes in the RS model in 
e + e~ and hadron hadron colliders at leading order in the gravitational coupling k. The leading 
order calculation requires the Feynman rules at order k 2 . This process displays a resonance 
structure due to the contribution from the s-channel KK mode exchange. For a large range of 
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Figure 4: The differential cross section for ^J~s = 2 TeV as a function of cos 9 for the production 
of the first KK mode in e + e~ collisions. Here 9 is the center of mass scattering angle and we 
have fixed c = 0.01, m = 150 GeV. 

the parameter space this resonance structure is visible in the two graviton production channel. 
The gravitons in the final state can be identified through their decay products which include 
a pair of jets, a lepton pair or a vector boson pair. We find that LHC can study a fairly 
large range of the parameter space of the RS model with this process. The process provides a 
nontrivial test of the RS model since its calculation requires the next to leading order Feynman 
rules. 
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Figure 5: The invariant mass W 2 dependence of the two graviton production cross section in 
proton proton collisions at = 14 TeV for several different values of the parameters (co, mo), 
where mo is given in GeV. Here W 2 = ts is the invariant mass of the two graviton final state. 
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Figure 8: The invariant mass W 2 dependence of the two graviton production cross section 
in proton anti-proton collisions at y/s = 2 TeV for several different values of the parameters 
(co,m ), where m is given in GeV. Here W 2 = ts is the invariant mass of the two graviton 
final state. 

Appendix A: Coupling to Matter Fields 
B.l Coupling to Scalar Bosons 

For a general complex scalar field $,we have 



C,(g) = g pa D p &D a <S>-ml&<f> 



tig 
5 2 £< 



\1V 



g=T) 



= 



g=T) 



5g^ u 5gP a 

where the gauge covariant derivative is defined as 

D^ = d, + igA;T a 
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Figure 9: The cos 8 dependence of the two graviton production cross section in proton anti- 
proton collisions at = 2 TeV. Here 9 is the center of mass scattering angle in the partonic 
subsystem. 

with g the coupling, the gauge fields and T a the Lie algebra generators. 

Substituting the expressions as given in Eq. (|20|) (j21j) and Eq. (j22j)in the lagrangian given 
in Eq. (fTTjl . we finally have the Feynman rule for the vertex involving two gravitons, and two 
scalars as follows: 

hfyh^QQ '■ llVwihfihp + h p k 2 p) + V^(k lp k 2u + k lv k 2p ) + f] vp {k XiJ k 2u + hahp)} 

+ [Vfipihahu + huha) + V^ihphcr + k la k 2p ) - Vpaihfihv + huhp)} 

- (VppVua + VvpVpv - VvvVpv)(h-k2 - ml) (24) 

This can be rearranged in the following simplified manner, 

^~\2 (Cpv,pa m Q 4" Cp,v, palx-qkik^fj (25) 

where k\, k 2 are the four-momenta of the scalars, 
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Figure 10: The rapidity dependence of the two graviton production cross section in proton 
anti-proton collisions at y^s = 2 TeV. 

and 

C^i/^palxrj — \j1fi\C pa, vr] f]cr\C ' pv, pr/ 1]p\Cpu, or\ f]v\C px), per 1]\rjCp Ut p a + (A <— > ?])] . 

Similarly, one can obtain the Feynman rules for the vertices involving two gravitons, two scalars 
and one gauge bosons(and two gauge bosons) as follows: 

h$h<ffA%<M> : pa Uh + k 2 yT a MN (26) 

7T 

h$hWA$A b x ** : -jpC m pa \ X v{T a , T b } MN (27) 



B.2 Couplings To Fermions 

Vierbein formalism is required in order to incorporate the fermion in the gravitation theory. 
The Fermion lagrangian is given by 

A/>(<?) = i)il P V P i) - m^ipip = -(ipiY'Dpip - VpipiYip) - m^ipip (28) 
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The covariant derivative V p i/j is given by, 



= (D p + -wfa ab )4> 



(29) 



The spin connection, w p b can be written as, 



Wf.ab = ~ (d p e bu - d v e hl j) e v a - - (d p e au - d u e ap ) e\ - (d p e ca - d a e cp ) (30) 

where e a p e b u r] ab = g pu , Y = e^ a , e%e bp , = rj ab and a ab = \[y^7b] ■ The conserved energy- 
momentum tensor of the fermion field is given by, 

T tv = -rj P u^il a D a ^ - m^ip) + l - U^ p D u ^ + ip^D^)] + ^q^^f^) 



(31) 



Now taking the functional derivative of C^{g) successively, one finds that 



5g^ u 5g(> a 



7 T 7 if;) 



9=11 



128 



+ V,ad a {^e ypaTl T i^) + Vuad Q (^e ppaTl T r^) + ^(V'WtVVO] (32) 



After substituting Eq. (|3T1) and in the expression of C^ n (Eq. (fT7}0 . we get the 

Feynman rule for the vertex involving two fermions and two gravitons as follows, 



htyh^W ■ tt^<M[?M7p(Pi + P2)u + l»(Pi + P2) P - 2vW +1/2- 2TJty))] 

+ [f]pp(la(Pl + V2)v + 7f(Pl + V2)a ~ + P2 ~ 2*71^,))] 

+[^p(Tm(Pi + Vi)<j + 7<r(Pi + - 2r/ M(7 (^/i +/ 2 - 2m^))] 
+ [^(7m(Pi + Ps)p + 7p(Pi + P2)^ - 2^ w (pi +1/2- 2m^))] 
-[^(7p(Pi + P2) CT + 7<r(Pi + P2) P - 2^(pi + y£ - 2m^))] 
-[^(Tm(pi + + T«/(Pi + P2V - 2^ (pi + p2 - 2m^))] 
+2C Mi ,, p CT [(^ +P2 - 2m,/,)] + -[?7 Aip e l , (T Ar7 T 7 5 (Pi ~P2) A 
+r]^pXrl T l 5 (Pi ~ P2Y + V^p P XrYl 5 (Pi - P2Y + Vupe^XrYl 5 (Pi ~ P2) X ]] (33) 

where 7 5 = 7°7 1 7 2 7 3 and Smn is defined in the flavour basis (M, N are the flavour indices) and 
Pi, P2 are the four momenta of the two fermions. 
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Similarly the Feynman rules for the vertices involving two gravitons, two fermions and one 
gauge particle is, 



A! 



'HP, \r) 



VppVXr,) 



+l v [{Vvp(Cpo, \v - VpvVxv) + ( u <-> A*)}] 

+Y' YlapiCpv, \ V - VpuVXv) + Vp.v{Cpa, Xr, ~ VpaVXr,)} 



(34) 



B.3 Coupling To Gauge Fields 

The lagrangian and conserved energy-momentum tensor for the gauge field are given by, 



C A (g) = ~F^F TV + \m A 2 A^A T - ±(d T A T + T; x A x ) 2 



(35) 



-F TT1 F TV - -m A 2 A T A T 
4 1 2 



F;F VT - m 2 A A^A v \ + - [d^d T A T A u + d v d T A T A, 



:Vpu 



d v d T A T A^ + -(d T A T : 



After successive functional differentiation, one obtain, 



s 2 c, 



5g^ u 5g 



pa 



9=11 



~ [F^F vp + F^F ua } s=v + [d p d^A u A a + d p d u A^A a ) 
+^ [d v dpA v A p + d (T d u A fl Ap\ - ^- Vpa [(d,A u + d u A fl )d v A n ] 
~^n P a l(d v d Pi A u A n + d^ApA")] 



1 

44' 



(d T d v A v A T + (d T A T ) 2 



(36) 



(37) 



The sets of Eq. fj3*Hj) . (|3l)|) and (|3Tj) together when substituted in Eq. (fTTjl . gives us the 
Feynman rule for the vertex involving two gravitons and two gauge fields as follows: 



hfttfPKK ■ J2 6ab [{v^(C P », r V (h.k 2 + m\) + D pVi TV (k h k 2 ) + r'Epu, r V ) + {p~<r) 
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+ 



Vu P (C all , r V (h-k 2 + m 2 A ) + D afJi: TJ) (k u k 2 ) + f 1 E a ^ TJ) ) + (p <-> <r) 



per, r?7 



per, r?7 y 



/l^, per 



(fci.fc 2 + m 2 A )r] TV - (k lrj k 2r - £ 1 k lr k 2v ) 

+ [(Vurilpvkipha + rjuqVrpki^p) - (r] UT 7] av k lp k 2p + 7] uv r] Ta k lp k 2p ) + (p <-> cr)] 

1 1 

-7 [VurVa V (k lp k lp + k 2p k 2p ) + (p <-> i/)] - - [ri VT ri pr) (k lp k la + fc 2M fc 2<T ) + (p <-> 

+^ [i]iyrVpa{kivk2p + k lp k 2v ) + (p <-> 1/)] + J [i] UT i] pa (k lrj k lfl + k 2p k 2r) ) + (p <-> 1/)] 
-| [C^, P a{ki v k lT + k 2r k 2v ) + (p <-> 1/)] - J [2C^ ^(^lr^)] 



(38) 



where fci and rc 2 are the four-momenta of the gauge fields A T and A^. 

Proceeding in the same way, one also finds the Feynman rules for the vertices involving two 
gravitons, three gauge bosons and two gravitons, four gauge bosons 



as follows, 



- -r2-/ abc [[{^C w , Au(fci - fa)* + VauC PPt \ S (k 3 - fci)^ + r] au C pp> vS (k 2 - k 3 ) x + (p <-> 

+ {(P O-)}] + {VtrvF m Xr)S (k u h, fc 3 ) + Aitf &2, fcs) + (P <-»• °~)} 

-VpaC pUy Xv (ki - k 2 ) s - i]paC pU) \ S {k 3 - h) v - i] pa C pU) vS {k 2 - k 3 ) x 
—VpisCpa, x n (ki — k 2 )s — i] pi/ C p(T: \s(k 3 — ki) v — i] pi/ C po - t v s(k 2 — k 3 ) x 
-f]po F pv, x v s(h, k 2 , k 3 ) - r] pu F p ^ Xv s(h, k 2 , k 3 ) 
+C pu , P a(vxv(ki - h)s + i] vS (k 2 - k 3 ) x + rj XS (k 3 - h) v ) 
-[{(fa/SIs* ~ k 3a r] pS )(r]priVuX - Vu V V P x) + Vus(Va V Vpxk 2p - rjpvVaxh^ + r} pr] r} pX k la - r} pX r} m k 2a ) 

+Vp8(VaXVu V k 2p - Va V Vuxh p + VpvVuxha - VlJ,XVu V k 2 a) + (p <-> cr)} - (v <-> P, P <-> o")]](39) 

Finally the Feynman rule for the vertices involving two gravitons and four gauge bosons 
are given by, 



-^■[{Vavir^f^Gpp, X V 5t 



reab recdf~i 
J J ^ftp, XSrjT 
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f ead r c G pp , Xv t5) + " ")} 



+ {Vu P {r C f ebd G^ X V 8r + r b r d G^ XStjt + r d f ebC G^ Xrir5 ) + (// <- u)} 

/ /eac rebds-t , /eab recdsi , /ead rebc/^i \ 

~Vpcr{j J ^pv, XifSr T/ J t-r^, Ao^r + J J <-»>!/, Ar;-n5j 

„_ / /eac rebdr~i , f eafe recd/^i , read rebc/^i \ 

—Vn^KJ J ^pcr, A»;or + J J <J><r, A<5r/r + J J <Jpcr, AryroJ 

"C^, P .{r C f eW {^mr - VXrVSr,) + rT^Vvr ~ VXrVSr,) + F^F^VtS ~ VXSVrr,)} 

-[{f eab f ecd (VvvVvxV P 8 - VpSVvxVvv) + f eca f ebd (v m VvxV8a - VpvVvSVvx) 

+ f ebC f eda (VpXVuvVSa - VpXVar,Vus)}%r + {f^ f^iVliSVuXVav ~ Vp5Vur,Vax) 
+f eCa f ebd (Vu8V m VaX - VmVvXVSa) + f^r^iVpXV^VSa - %XVa V Vus)}Vpr 
+ {f eab f £Cd (V8aVpvVpX ~ VSvVpXVpv) + f eCa f ebd {VpSVpXVa V - VpSVpXVar,) 

+ f ebC f eda (Vp8VpvV*X - Vp8Vpr)V<rx)}VvT 
+ {f eab f eCd (Vu8VpXVpr, - VuSVpvVpx) + f eCa f ebd (Vp8VpXV,r, ~ VpSVpXVu V ) 

+f ebc f eda (VpSVpvVuX - VpSVpvVux)}Vt8 + {p^ v}} (40) 

where, 

p<j(k u k 2 ) = Vpvhahp - [Vp<rhuk2p + Vpphahv ~ Vptrfa^v + (/i <-> f)] 
^Tif, p<7(ki,k 2 ) = V^(klpkla + &2p^2a + hpha) ~ [Viyahphp + Vvpfopfav + ^ ^)] 

,KrA(fci,fc2,fcs) = VppVaxih - h) u + rj^ripxih - k 1 ) v + ij^VpAh ~ k 2)v + (fJ> <-> 

C^, pCT A5 = V^iVpaVXS - VpSVctx) + [VppVuSVXa + VpXVuaVpS - VppVuaVXS - V pX^SV pa] 

Appendix B: Triple graviton vertex 

To find the feynman rules for three graviton vertex we expand the pure gravity Lagrangian 
upto next order in weak field approximation. In this approximation 

L g = \^fgR = L Q + kL ± + ... 

= 1 (d^ p d,h vp - d»Kd,h p p - d^d p h pv + 2d,h^d u hp) 



4 

K 

4 

K 

4 



2h^d^h f , a d a h uP + h" u d a h^d a h tiu - 2h^d a h lx pd a hP] + .... (41) 
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From the expression given in L 1; feynman rule for three graviton vertex reads 

g [(^2m^3i/ + k 2v k 3p )F p xka + {k\ p k 3(J + k la k 3p ) F^xki* + (k ix k 2k + k lk k 2X ) F ppau ] — 
g [(k 2p k 3u + k 2l/ k 3p )r] Xk r] pa + V^Vxkihpha + k la k 3p )} - 

- [r)pvVpa(ki\k 2k + k lk k 2X )} + 

o 

g [Fp P \vk 2 kk 3o - + F ppkv k 2 xk 3a + F pa xuk 2k k 3p + F pakv k 2 \k 3p \ + 

■g" [-^ipcrA^lfc^ + Fp,pakk\\k 3v + F upa \ki k k 3p + -Fl/po-fc^U^u] + 

— [F p x kp ki a k 2u + F p xkaki p k 2u + F u \ kp ki a k 2u + F v \ ka k\ p k 2p \ + 

— [G \iv P \ckk\-k, 2 + G pap \ vk k 2 .k 3 + Gp^Apfco-^i-^] — 

- [/ Ppff A^il 2 + ^2-A;3 + fc 3 -A;i)] (42) 
where F, G and I are defined by the following relations 

^r/ p(T + 77^77^ (43) 

G pup \ak = VpvFpXka + f]paF p Xku (44) 
I p,vpa\k T]ppF u Xk<T ^jpaFuXkp ~\~ ^]vpF p Xka ~\~ ^]vaF p Xkp 

(45) 



Appendix C: Cross section for two graviton production through fermion 
anti-fermion annihilation 



In this Appendix we give the detailed form of the cross section for two graviton production 
by fermion anti-fermion annihilation. The differential cross section for this process is given by. 



da 



1 4M 2 

■ r = V 1 ^(Tii+T 22 + T33 + T 12 + T 2 3 + T 1 3+T44 + T 1 4 + T 2 4 + T 3 4) (46) 

d(cosO) 32ns V s 



where the terms inside bracket are given by, 

7r2 c 4 /-2(M n 2 - tf (9M n 8 - 4t 3 u + 12M n 2 t 2 (t + u)- M n H (20* + 9m)) N 



Tn = 



9M 8 t 2 



(47) 



20 




Figure 11: A Typical Triple Graviton Vertex. 

7t2 c 4 /-2(M„ 2 - u) 4 (9M n 8 - Atu 3 + 12M n V (t + u) - M n \ (9t + 20m) 
T22 = ~^4\ 9M n V 



+ 



12 



7r2 c 4 /2 (T092M n 12 - 1230M n iU (t + u)J + 16t(t - u) u(* + u) 

tt 2 c 4 / 7M n 6 (t + u) (25t 2 + 52tw + 25m 2 ) \ 

tt 2 c 4 /4M n 8 (113t 2 + 299tw + 113m 2 ) \ 

tt 2 c 4 (M 2 (t + u) (15t 4 + 62t 3 M - 34t 2 w 2 + 62tw 3 + 15w 4 )' 
~K V 9M^ , 

ir 2 4 ( M n A (56t 4 + 370t 3 M + 84tV + 370tw 3 + 56m 4 ) \ 
~< V 9M? J 

^-2 C 4 / (2M„ 2 - t - w) (88M n 14 - 62M„ 12 (t + u) - 4t 3 w 3 (t + u 

9M^tu 

/ (2M„ 2 - t - m) (lOM n 10 (7t 2 + 19tu + 7m 2 )) ' 
^ 9M^tu 
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7r2 c 4 / {2M 2 -t-u) 2M n 2 t 2 u 2 (9t 2 + Utu + 9u 2 ) \ 
' 9MHu 



4 



7l 2 C 



m} 



.4 / (2M n 2 -t-u) (M n Hu (t + u) (26t 2 + 47tw + 26m ; 



9MHu 



+ 



•2 C 4 / (2M r 2 -t-u) (M n 8 (t + u) (62t 2 + 41tu + 62m 2 )) ' 



71 



9MHu 



+ 



7^4 / (2M„ 2 -t-u) (M n 6 (— 12t 4 + 64t 3 M + 202t 2 M 2 + 64tu 3 - \2u 
' 9M n 8 tu 



4 

ml 



ml 



ml 



4 

m 



2^4 



— 99M n 14 + 2t 4 (t — u)u (t + u) + M n 12 (96t + 63m)) \ 



9M*t 



+ 



M n 10 (5t 2 + 70tu - 9m 2 ) - M n H (3t 2 + 92tu + 55m 2 )) 



9M*t 



M n H 3 (3t 3 + 12t 2 M - htu 2 + 2m 3 ) + M n H (-18t 3 - 40t 2 M + tu 2 + 6m 3 ) 



7T C, 







9M n H 

M n H 2 (13t 3 + 36t 2 M + 31tM 2 + 7m 3 ^ 



7T^ 
4 



2^4 



9M„ 8 t 



-99M n i4 -2t(t- u) u (t + u) + M n u (63t + 96m 



7T C, 







9M n m 

M n 8 u (55t 2 + 92tM + 3m 2 ) + M n 10 (-9t 2 + 70tM + 5m 2 



9M„ 8 m 



+ 



7r 2 c 4 /8 



M n 6 M (6t 3 + t 2 M - 40tM 2 - 18m 3 ) - M n V (2t 3 - 5t 2 M + 12tM 2 + 3m 3 )) 



2^4 



7T C, 







9M n u 

M n 4 M 2 (7t 3 + 31t 2 M + 36tM 2 + 13m 3 



ml 



9M*u 



X 2 cl / 8688M W 16 - 9720M TO 14 (t + u) + 9t(t - m) 2 m(£ + m) 4 \ 
44 " 512^1 V 9M? J + 
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X 2 s Cq (20M n M (t + u) (85t 2 + 52£m + 85m 2 ) \ 



512m 8 , V 9M„ 8 

\ 2 s ct ( M n 12 (428t 2 + 5048tM + 428m 2 ) \ _ 
512mf V 9M n 8 J ~ 

A 2 c 4 ( YlM 2 (t + uf (t 4 + 8t 3 M - lOtV + 8tM 3 + m 4 ) \ 
512m 8 , V 9M„ 8 / + 

A 2 c 4 / 4M n 8 (17t 4 + 403t 3 M + 24tV + 403tM 3 + 17m 4 ) \ _ 
512m 8 , V 9M n 8 / 

X 2 4 ( AM n & (t + m) (117t 4 + 482t 3 M + 410t 2 M 2 + 482tM 3 + 117m 4 ) 
512mf V 9M n 8 
A 2 c 4 /M n 4 (t + uf (135t 4 + 604t 3 M + 58tV + 604tM 3 + 135m 4 ) \ 
512m 8 , V 9M„ 8 J 



ttc 4 A s ^864M n 16 + 3t 4 (t-M)w(t + M) 2 \ 



+ 



4mg V 9M r °t 
7t4\ s (2M n H 4 (3t -u)(t + u) (t + 5m) 
Ami V 9M„ 8 t 
ttc 4 A s / 12M W 14 (80t + 63m) + M w 12 (-79t 2 - 404tM + 207m 2 ) 
Ami V 9M n 8 t 
7TCqA 5 / M n 8 t (39t 3 - 49t 2 M - 323tM 2 - 171m 3 ) \ _ 
Ami V 9M n 8 t J ~ 

ttc 4 A, / 2M n 10 (153t 3 + 593t 2 M + 331tM 2 - 9m 3 ) \ 
Imf" V 9M n 8 t J + 

7TCqA s / 2M n 6 t (67t 4 + 201t 3 M + 175tV + 37tM 3 - 6m 4 ) \ 
Ami V 9M n H J ~~ 

71CqX s (M n H 2 (51t 4 + 186t 3 M + 143tV + 94tM 3 + 18m 4 ) \ 
Ami V 9M„ 8 t J 



vrc 4 A s / (864M ra 16 - 3t (t-M)M 4 (t + M) 2 )' 



Ami I 9M„ 8 m 

vrc 4 A s (2M 2 (t - 3m) m 4 (t + u) (5f + u) 
Ami V 9M n 8 M 
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+ 



+ 



ttCq\ s ( 12M W 14 (63t + 80m) + M n 12 (207t 2 - AOAtu - 79m 2 ) 

4mg V 9M n 8 u 
ttcqX s (M n 8 u (-17K 3 - 323t 2 M - 49tM 2 + 39m 3 ) \ 

4mg V 9M n 8 w / 

7TCgA s /M n 10 (-18t 3 + 662t 2 M + 1186tw 2 + 306m 3 ) \ 

Ami V 9M„ 8 « / " 

%Cq\ s (M n \ 2 (18t 4 + 94t 3 M + 143t 2 M 2 + 186tM 3 + 51m 4 ) \ 

Ami V 9M„ 8 m ) 

ttc 4 A s /2M n 6 M (-6t 4 + 37t 3 M + 175t 2 M 2 + 201£m 3 + 67m 4 )) 

Ami \ 9M n 8 M 



+ 

) (55) 



vrc 4 A s / 16 (-588M n 14 + 630M n 12 (t + u y 
Tu ~ 32^§ I 9lC 



7TCqA s / 3t(t — u) 2 u(t + m) 



32m[j V 9M n 8 J 

ttc^\ s ( M n 8 (t + m) (5t 2 - 124tM + 5m 2 ) ' 



32ml \ 9M n 8 J 

7rc 4 A s /2M n 10 (43t 2 + 244tM + 43m 2 ) \ 

32ml \ 9M n 8 J ' 

%4\ s (M 2 (t + uf (3t 4 + 22t 3 M - 26t 2 M 2 + 22tM 3 + 3m 4 )' 

32ml \ 9M n 8 

ttc 4 A s /4M„ 6 (8t 4 + 61t 3 M + 30t 2 M 2 + 6Um 3 + 8m 4 

32ml \ 9M n 8 / 

ttc 4 A s /M n 4 (t + m) (21t 4 + 92t 3 M + 38t 2 M 2 + 92£m 3 + 21m 4 ) 

32///!; V !).!/„ 



+ 

(56) 



t \ 

Here s, t and u are the Mandelstam variables. 
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